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We investigate quantum correlations in the ground state of noninteracting Fermi gases of N parti- 
cles trapped by an external space-dependent harmonic potential, in any dimension. For this purpose, 
we compute one-particle correlations, particle fluctuations and bipartite entanglement entropies of 
extended space regions, and study their large- N scaling behaviors. The half-space von Neumann 
entanglement entropy is computed for any dimension, obtaining Sns ~ ciN^ d '~ 1 '' d IniV, analogously 
to homogenous systems, with ci = 1/6, l/(6v2), l/(6v6) in one, two and three dimensions respec- 
tively. We show that the asymptotic large-iV relation Sa ~ 7t 2 Va/3, between the von Neumann 
entanglement entropy Sa and particle variance Va of an extended space region A, holds for any 
subsystem A and in any dimension, analogously to homogeneous noninteracting Fermi gases. 



PACS numbers: 03.65.Ud. 05.30.Fk, 03.67.Mn 



I. INTRODUCTION 

The recent developments in the experiments of dilute 
atom gases, namely the achievement of Bose-Einstein 
condensation in dilute atomic vapors of 87 Rb and 23 Na (J 
2] and the great progress in the experimental manipula- 
tion of cold atoms in optical lattices 0] have provided a 
great opportunity to investigate the interplay between 
quantum and statistical behaviors in many-body sys- 
tems. A common feature of these experiments is the 
presence of a confining harmonic potential which traps 
the particles within a limited spatial region. The capa- 
bility of varying the confining potential, which may also 
depend on the spatial directions, allows also to vary the 
effective spatial geometry of the particle systems, includ- 
ing quasi-lD geometries, see, e.g., Refs. 

In this paper we investigate the quantum correlations 
arising within the ground state of noninteracting Fermi 
gases trapped by an external space-dependent harmonic 
potential. Quantum correlations can be characterized by 
the expectation values of the products of local operators, 
such as the particle density and one-particle operators, 
or by their integral over a space region A, such as the 
particle-number fluctuations within A. Quantum corre- 
lations are also characterized by the fundamental phe- 
nomenon of entanglement, which gives rise to nontrivial 
connections between different parts of extended quantum 
systems flol - fl2l |. A measure of entanglement is achieved 
by computing von Neumann (vN) or Renyi entanglement 
entropies of the reduced density matrix of a subsystem. 
One-particle correlations and bipartite entanglement en- 
tropies provide important and complemetary information 
of the quantum behavior of many-body systems, because 
they probe different features of the quantum dynamics. 

We consider Fermi gases of N particles confined by 
harmonic traps of arbitrary dimension, and study the 
large- N scaling behavior of the above-mentioned observ- 
ables to characterize the quantum correlations of the 
ground state. We determine the asymptotic behaviors of 
the half-space entanglement entropies in any dimension, 
which turn out to increase as N^-V/dlnN, analogously 



to homogenous systems [Hj]. We study the relation be- 
tween particle fluctuations and entanglement entropies 
of extended space regions. This is motivated by recent 
proposals of considering the particle fluctuations as ef- 
fective probes of many-body entanglement at zero tem- 
perature fl4rfl9j , which are more easily accessible exper- 
imentally. In homogeneous finite- volume systems of non- 
interacting fermions, of any dimension, the vN entangle- 
ment entropy Sa of an extended subsystem A turns out 
to be closely related to the particle variance Va within 
A. Indeed, asymptotically for a large number of parti- 
cles N, Sa/Va ~ 7r 2 /3 for any subsystem A and in any 
dimension 19], with 0(1/ IniV) corrections. We show 
that this asymptotic behavior also holds in the presence 
of a space-dependent harmonic potential, such as the one 
which characterizes recent experiments with cold atoms. 
For this purpose we present several analytical results in 
the large-iV limit, and numerical (practically exact) re- 
sults at fixed N by computations from the ground-state 
wave function. 

The paper is organized as follows. Sec. [TT] reports some 
general expressions for the ground-state many-body wave 
function of free fermion gases in a harmonic trap, and 
define the observables that we consider. Sec. IIIII focuses 
on one-dimensional (ID) systems. Systems in higher di- 
mensions are considered in Sec. IIVI Finally, in Sec. IVl we 
summarize our main results and draw our conclusions. 



II. GROUND STATE AND OBSERVABLES OF 
TRAPPED FREE FERMION GASES 

A. The ground state in a harmonic trap 

We consider a gas of N noninteracting spinless 
fermionic particles of mass m confined within a limited 
space region by an external potential. In the following 
we set H = 1 and m = 1. The ground-state wave function 
is 



*(xi, ...,xjv) 



1 



det [■;/>; (xj)], 



(1) 
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where ipi are the lowest ./V eigensolutions of the one- 
particle Schrodinger equation 



HrP^EiiPi, H = ^- + V{x). 



(2) 



A generic power-law rotational-invariant potential such 
as 



V(x) 



1 /x 



2\ P/2 



P V i 2 



(3) 



where ^ is the trap size, gives rise to a trap length scale 
£ which behaves as a nontrivial power of Z, 



p 



p 



(4) 



where 9 is the trap exponent [20j , which does not depend 
on the spatial dimension in free fermion gases. The power 
p = 2 describes the harmonic trap, where £ is the so- 
called oscillator length. In the limit p — > oo the system 
becomes equivalent to a Fermi gas confined by a hard- 
wall spherical trap of radius I. 

The one-particle energy spectrum in harmonic traps is 
discrete. The eigensolutions can be written as a product 
of eigenfunctions of corresponding ID Schrodinger prob- 
lems, i.e. 



'/ , Tii,ii2,...,nj( x ) — J_ J_ fin, (^i); 



•Em, 



8=1 
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where the rii label the eigenfunctions along the d direc- 
tions, which are 



i(*) = r 1/a - 

TT- 1 - 



H n -!{X) 



7r l/4 2 (n-l)/2( n _ 1)11/2 



-X 2 /2 



r 2 (n- 1/2), n=l,2,.. 
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where X = x/£, and H n {x) are the Hermite polynomi- 
als. Note however that, although the spatial dependence 
of the one-particle eigenfunctions is decoupled along the 
various directions, fermion gases in different dimensions 
present notable differences due to the nontrivial filling of 
the lowest N states which provides the ground state of 
the iV-particle system. 

The above one-particle eigensolutions allow us to re- 
construct the corresponding fermion-gas ground state 
(HJ), and study its general properties by computing parti- 
cle correlations and bipartite entanglement entropies. In 
the following we also set 1 = 1, thus 



£ = Z e = l, X = x/£ 



(9) 



The dependence on fi, m and I = cu^ 1 of the quantities 
considered can be easily reconstructed by a dimensional 
analysis. 



Observables 



1. One-particle and density correlations 



The one-particle correlation function reads 

JV 

C(x,x) = (ct(x)c(y)> = X>(x)*lMy) (10) 

i=l 

where c(x) is the fermionic annihilation operator. The 
particle density and the connected density-density corre- 
lation are respectively given by 

N 

p(x)Mn(x))=C(x,x)=]ThMx)| 2 , (11) 

i=i 

G„(x,y) = (n(x)n(y)) c = 

= -|C(x, y )| 2 +£(x- y )C(x,v) (12) 

where n(x) = c(x)^c(x) is the particle-density operator, 
and we used the Wick theorem to write G n in terms of 
the two-point function C . 



2. Particle fluctuations and entanglement entropies 

Other important measures of the quantum correlations 
are related to extended spatial regions, such as the dis- 
tribution of the particle number and the entanglement 
with the rest of the system. 

The expectation value and connected correlators 

„ rn m 

N A = (N A ) 7 (N%) c = / n^dJ^x^c, (13) 

JA t=l 1=1 

of the particle-number operator of an extended region A, 



N A = An(x), 



(14) 



characterize the particle distribution within A. For this 
purpose, it is convenient to introduce the cumulants of 
the particle distribution, which can be defined through a 



generator function as [21 



Vi m) = (-id x ) m ]n{e iX * A }\x=o. 



In particular, the particle variance reads 

(2) ,„ r 2\ _ , AT 2\ 
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(16) 



(the superscript m — 2 will be understood in the case of 
the particle variance). A measure of the entanglement 
of the extended region A with the rest of the system is 
provided by the Renyi entanglement entropies, defined as 



5' 



(a) 



1 — a 



lnTr P 2 



(17) 
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where pa is the reduced density matrix pa of the subsys- 
tem A. For a — > 1, we recover the vN definition 



S A EE S« EE 



-Tr In pa 



(18) 



(the superscript a = 1 will be understood in the case of 
the vN entanglement entropy). 

In noninteracting Fermi gases the particle cumulants 
and the entanglement entropies of a subsystem A can 
be related to the two-point function C(x,y) restricted 
within A, which we denote by Ca(x,u). The particle 
number and cumulants within A can be derived using 
the relations (see e.g. Ref. [l8l |) 



Na = Tr Ca , 



(m) 



{-id z ) m G(z,C A )\ z = , 



£(z,X) = Trln [l + (e 4Z - l) X] 



(19) 
(20) 
(21) 



The vN and Renyi entanglement entropies can be eval- 
uated from the eigevalues of Ca(x,v) (see Refs. [22I I23I] 
for applications to lattice systems). 

The computation of particle cumulants and entangle- 
ment entropies in Fermi gases of N particles is much 
simplified by introducing the N x N overlap matrix 
A ii, 



,,„„ ■ , d d zi>* n {z)i}) m {z), n,m = l,...,N, (22) 

I A 

where the integration is over the spatial region A, and 
involves the lowest N energy levels. The overlap matrix 
A and the restricted two-point function <Ca satisfy 



Tr(rt=TrA fc V k e N, 



(23) 



which implies that the particle cumulants and the entan- 
glement entropies can be computed form the eigenvalues 
of the N x N overlap matrix A. The eigenvalues a, of A 
are real and limited, a, € (0, 1). 

The particle number and cumulants can be computed 
using the relations 19] Na — TrA and 



v[ m) = (-id z ) m g(z,A)\ z = . 



In particular, 



V A = TrA(l — A), 

Vf ] =Tr[A-3A 2 + 2A 3 ], 

,(4) 



(24) 

(25) 
(26) 



V\ > = Tr[A— 7A 2 + 12A S - 6A 4 ] , (27) 

etc.... The vN and Renyi entanglement entropies are ob- 
tained by [H [H 



N 



&A ^ — ^ ] s a(°"n). 



where a n are the eigenvalues of A, and 
1 



s Q (A) 



1 — a 



In [A Q + (1 — A)"] 



(28) 



(29) 



and, in particular, 



Si(A) = — A In A - (1 - A) ln(l - A) 



(30) 



for the vN entropy. We also mention that, while TrA 
gives the average particle number Na within A, detA is 
the probability to find all particles within A. 

We consider two different partitions of the space: 
(i) The subsystem B is separated from the rest by 
a hyperplane at a distance x from the center of the 
trap. The corresponding x-dependent particle cumulants 
and entanglement entropies are denoted by Vjf 1 ^ (x) and 

Sg(x). In one dimension, the subsystem B is given by 
the infinite interval B = [—00, a;]. The half-space quanti- 
ties are defined as 



Vkf EE Vi" 



(0), s£j = s£\o) 



We also define 



(«), 



(31) 



(32) 



(ii) The subsystem S is a region containing the cen- 
ter of trap, and enclosed by two parallel hyperplanes at 
distance x from the center. The corresponding entangle- 
ment entropies are denoted by V$ (x) and s[?\x). In 
one dimension, the subsystem S is given by the symmet- 
ric interval S — [— x, x] (where x — corresponds to the 
center of the trap). 



III. FERMI GASES IN ID TRAPS 

In this section we consider ID noninteracting spinless 
fermion gases of N particles confined by a power-law ex- 
ternal potential, in particular by a harmonic potential. 
This model has a wider application, because ID Bose 
gases in the limit of strong short-ranged repulsive in- 
teractions can be mapped into a spinless fermion gas. 
The basic model to describe the many-body features of 
a boson gas confined to an effective ID geometry is the 
Lieb-Liniger model with an effective two-particle repul- 
sive contact interaction 



27 



/v 



2m 



V( Xi 



where N is the number of particles and V(x) is the con- 
fining potential. The limit of infinitely strong repulsive 
interactions corresponds to a ID gas of impenetrable 
bosons [28[, the Tonks-Girardeau gas. ID Bose gases 
with repulsive two-particle short-ranged interactions be- 
come more and more nonideal with decreasing the par- 
ticle density, acquiring fermion-like properties, so that 
the ID gas of impenetrable bosons is expected to pro- 
vide an effective description of the low-density regime of 
confined ID bosonic gases [29]. Therefore, due to the 
mapping between ID gases of impenetrable bosons and 



4 



spinless fermions, some correlations in free fermion gases 
are identical to those of the hard-core boson gases, such 
as those related to the particle density, particle fluctu- 
ations of extended regions, and bipartite entanglement 
entropies of connected parts. This correspondence holds 
also in the presence of an external space-dependent po- 
tential. 



A. The particle correlators 

In ID noninteracting Fermi systems with N particles 
in a harmonic trap, the two-point correlation function 
(fT0|) can be written as 

ru x N 1 / 2 (j) N+1 (x)(j) N (y) - (f> N (x)(j) N+1 (y) 

c \ x >y) = —/K > ( 33 ) 



0.0 



V2 



x-y 



where we used the Christoffel-Darboux relation for oth- 
ornormal polynomials. The particle density p(x) = 
C(x,x), 



p(x) 



TV 1 / 2 
V2 



[(j)' N+1 {x)<p N (x) - (f>' N (x)(f> N+ i(x)] , (34) 



shows a peculiar behavior characterized by N local max- 
ima, which get suppressed by powers of 1/N with increas- 
ing N. 

Since the particle density and the density correlator of 
free fermion gases are equal to those of boson gases in the 
hard-core limit, the results already obtained for system of 
impenetrable bosons in a trapping potential apply also to 
trapped fermion gases. The large-TV asymptotic expan- 
sion is known to 0(1/N) [3(1 HH. The leading behavior 
is given by 



p(x)=N 1 / 2 [R p (0+O(l/N)], t^x/N 1 / 2 , (35) 



with 



(36) 



and R p (C) = for C > Cc = V2. 

The space dependence of the connected correlation 
function of the particle density operator n x = c(xy c{x) 
presents a different large- N scaling behavior, character- 
ized by different power laws 32J, i.e., 



G{x,y) nNR G {N l > 2 x,N l ' 2 y), 



(37) 



for x 7^ y, as shown by Fig. [T] Note that the asymptotic 
regime of G n is not approached uniformly when x — > y, 
because G(x,x) = p(x) — p(x) 2 , cf. Eq. (jl2]l . 

The large- N scaling function of the particle density in 
a harmonic trap vanishes at £ c = v2 where R P {C, C ) = 0. 
Around this point the particle correlations (32|, and also 
the entanglement entropies, show a different large- AT scal- 
ing behaviors characterized by other power laws [33| . In- 
deed, around the point x c = jV 1 / 2 ^, the particle density 
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FIG. 1: (Color online) Large- TV scaling of G n (x,y) for x = 
y/2, x = 0, and x — —y for ID fermion gases in a harmonic 
trap. We plot N^ 1 G„(x, y) vs N 1/2 y. 



and its correlation behave as 

p{x)^NV & g p [N l l\x-x c )}, x e =NV 2 C c , (38) 
G n {x c ,x) w N^g^N^ix - x c )]. (39) 

The scaling function g p (z) can be obtained from related 
computations within the Gaussian unitary ensembles of 
random matrices [3l|, HH : 

g p (z) = UmN^N-^plN^iCc + N-^z)} 

= 2 1 / 2 |Ai'(2 1 / 2 z)| 2 -2z|Ai(2 1 / 2 z)| 2 . (40) 

B. Spatial entanglement 

1. Half -space entanglement entropy 

The asymptotic large-TV behavior of the half-space 
(A = [—oo,0] where x = is the center of the trap) 
vN and Renyi entanglement entropies can be inferred 
by exploiting known results for the ID hard-core Bosc- 
Hubbard model in the presence of an external power-law 
potential V(x) = (x/l) p and a chemical potential, which 
is equivalent to a lattice free-fermion model. The deriva- 
tion is outlined in App. [5] We obtain 
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°ASY — u » 



o(N°), 

m/v + ini^ + yo 



where 



G a — 



12 



(41) 
(42) 

(43) 



p is the power-law of the potential, and y a is given in 
Eq. (|A12|) . Note that the leading logarithmic term, and 
in particular its coefficient, is independent of the trapping 
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FIG. 2: (Color online) The half-space vN entanglement en- 
tropy. The full line shows the large-iV asymptotic behavior 
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eigenvalues and then obtain the entanglement entropies 
through Eq. ([25]). Figs. H and [3] show data for the vN 
and a = 2 Renyi entropies. They are fully consistent 
with the asymptotic behavior (|4"2"1) . In particular, the 
large- N behavior of the vN entropy turns out to behave 
as 



'HS 



Sasy + ( — 1 



N 



C2 

iV 2 



+ (-1) 



jV c 3 

iV 3 



(45) 



with ci ~ —1/8 (with a precision better than 10 6 , see 
Fig. El c 2 « 0.009893 and c 3 0.066 (the uncertainty 
should be on the last figures). The data of the a = 2 
Renyi entropy, shown in Fig. 02 fits the Ansatz 



? (2) 
'hs 



? (2) 

'asy 



+ (-1) 



AT 



A/1/2 



62 

N 



(46) 



with 61 « -0.2387 and b 2 « 0.0146. Note that the above 
numerical results show that the corrections to the large- 
N asymptotic behavior are 0{N~ x / a ) in Eq. (14"T]) , anal- 
ogously to homogeneous systems, cf. Eq. (|44l) . 

Finally, we mention that the effects of a power-law 
trapping potential on the scaling behavior of the entan- 
glement at the quantum critical point of ID lattice mod- 
els were investigated in Refs. [33l 13614381 ]. In particular, 
the ID hard-core Bose-Hubbard model, which is equiv- 
alent to a free fermion lattice model, was considered in 
the superfluid phase at half filling [33[. As shown by the 
arguments reported in App. [3] used to derive Eq. (|42l) . 
these results are somehow related with the large- N scal- 
ing of ID Fermi gases investigated in this section, in 
particular when the chemical potential is driven toward 
the supcrnuid-to-empty transition. However the large- N 
scaling behavior of ID continuum Fermi gases presents 
distinct features, as pointed out in Refs. [2J, [26( in the 
case of homogenous systems. 



FIG. 3: (Color online) We plot (-1) 1+jv (Shs - Sasy) vs TV" 1 
(bottom) and (-1) 1+JV (S^ - Sfly) vs N~ 1/2 (top). In the 
bottom figure the dashed line shows the slope 1/8. 



potential, and it is equal to that of homogeneous systems 
with open boundary conditions [24], |26| , which is deter- 
mined by the corresponding conformal field theory 35 1 
with central charge c = 1 . The asymptotic behavior (|4"2")l 
in the limit p — > 00 reproduces the results for homoge- 
neous systems with open boundary conditions [24l l26j | 



S 



(a) 
HS 



lniV + ln4 + y Q + 0(N- 1/a ) 



(44) 



because in the limit p — > 00 the system becomes equiv- 
alent to a Fermi gas confined by a ID hard-wall trap of 
size L = 21. 

In order to check the convergence to this asymptotic 
behavior, we numerically compute the half-space entan- 
glement entropies S^g of N particles in the presence 
of a harmonic trap. We use the method based on the 
overlap matrix (|22[) . i.e. we numerically compute its 



2. Finite intervals around the center of the 



We now consider a symmetric interval S — [— x, x] 
around the center of the trap. By integrating the large- N 
particle density (|35j) within the interval S — [— x,x], we 
obtain the average number Ng(%) of particles within S 
in the large- TV limit, 



N s (x) 
N 



C\/2 - C 2 + 2arcsin(C/\ / 2) 



0(1/N) 

(47) 

Fig. 0] shows results for the particle 
100. They show 



where ( = x/N 1 / 2 . 
number within S at fixed N, up to N 
that the large- iV limit (I4T1) is rapidly approached by the 
data. 

Results for the vN and a — 2 entanglement entropies 
up to N = 180 are shown in Figs. [5] With increasing 

N, the subtracted data of — 2C a In N appear to ap- 
proach a function of C = x/N 1 / 2 . Therefore, we infer the 
large- N scaling behavior 



4 a) (x) w 2C a In N + y a + In 4 + f { s a) (Q 



(48) 
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FIG. 4: (Color online) The particle number within the inter- 
val S — [—x,x] around the center of the trap, for some values 
of TV up to TV = 100, versus ( = x/N 1/2 . The full line (hardly 
visible among the data symbols) shows the large- TV limit (|47[) 
of the ratio N s (x)/N. 

The scaling functions f s a ^ (C) are expected to be singular 
at £ = corresponding to a vanishing interval, and at £ = 
•\/2, which corresponds to the point where the particle 
density vanishes in the large- TV limit, cf. Eq. (|35l ). Note 
that the space dependence scales analogously to that of 
the particle density, cf . Eq. (pJ5j) , while it differs from that 
of the connected density-density correlation, cf Eq. (I57l) . 

The large- TV convergence is rapid at least up to £ m 1 , 
but also the data for 1 < ( < v2 appear to approach 
a unique curve, although more slowly. Moreover, the 
behavior of the data with increasing TV suggests that the 
large- TV scaling functions /£ (C) are independent of a. 
Actually, they turn out to be well approximated by the 
simple function 



fi a) (0 



f a (C) = lnsin«/\/2) + ln(2/7r), (49) 

as shown in Fig. [5] for the a = 1 vN and a = 2 Renyi 
entropies. In Fig. [5] we show the differences between the 
vN entropy Ss(x) and the asymptotic behavior (|4"5)> with 
fs given by Eq. (|49l) . at its maximum ( = 2 -1 / 2 and at 
( = 2~ 3 / 2 . For example at ( = 2 -1 / 2 the data show 
deviations smaller than 0.01 for TV > 100, suggesting 
that the deviation in the large-TV limit should be less 
than 0.01. Smaller deviations are observed at C, = 2 -3 / 2 , 
see Fig. [SJ A more precise large- TV extrapolation is made 
difficult by the presence of oscillations, whose structure 
is not clear, see Fig. 

It is worth comparing the above results with the be- 
havior of analogous quantities in homogeneous Fermi gas 
within hard walls, whose entanglement entropies of the 
interval S — [—x, x] around the center of the hard-wall 
trap of size L = 21 = 2 are given by 26] 



Sg i x ) ~ 2C Q 



lnTV + lnsin(7ra;) + y a + ln2 + 0(TV~ 1/q ) 

(50) 




FIG. 5: (Color online) The a = 1 vN (bottom) and a = 2 
Renyi (top) entanglement entropies of the interval S = [—a;, x] 
vesus C s x/N 1/2 . We plot (x)/(2C a ) - (lnN + y a + ln4) 
for some values of TV up to TV = 180. The full lines show the 
function (35) . 
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FIG. 7: (Color online) S^ (x)/C a versus :r/iV 1/2 for the a = 
1 vN (bottom) and a = 2 (top) entropies, and several values 
of jV. The two sets of data appear to approach the same 
large- N limit. 



Finally, Fig. [7] shows results for the quantity 5£(x) 
Sg(x) — S^ s , i.e. the difference between the entanglement 



entropies of the intervals [—00, x] and [— 00, 0]. 
show the large- N scaling behavior 



c Q /i Q) (C), 



C^x/N 1 / 2 . 



They 



(51) 



They also suggest that = 



is indepen- 



dent of a, although the convergence of the a — 2 Renyi 
entropy is slower than that of the vN entropy, apparently 
0{N- 1 / 2 ) against 0{N^ 1 ). 



C. Particle fluctuations in extended spatial 
subsystems 

Some results for the half-space particle variance and 
quartic cumulant are shown in Fig. [8] They are charac- 
terized by large odd-even oscillations in the number of 
particles. An educated guess for the asymptotic large- N 
behavior of the half-space particle variance in a generic 
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FIG. 9: (Color online) Check of the asymptotic behavior for 
the particle variance and quartic cumulant of half space. The 
subscript "av" indicates the average over the last two data to 
suppress the large oscillations. The dotted and dashed lines 
show the TV — > 00 limit expected for Vhs — Vasy ~ and 
v hs ~ -0.009255. 



external potential V(x) oc (x/l) p is 
Vhs « Vasy + o(N°), 



Vasy = 



1 

2^2 



In N + In 



4( P + 2) 
P 



1 



IE 



• (52) 



This asymptotic behavior is somehow derived by analogy 
with the asymptotic behavior of the half-space Renyi en- 
tanglement entropies, taking also into account the known 
asymptotic behavior of the particle variance in hard-wall 
traps [H, 



Vks 



1 

2^ 



[lniV + ln4+l+7, B + 0(jV- 1 )] , (53) 



which must be recovered in p — s- 00 limit. 
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Concerning the other cumulants, we expect that the 
leading term is the same as that of homogenous systems 
within hard walls, like the leading terms of the entangle- 
ment entropies and particle variance. Thus 



Vr 



(2i) 



HS 



V2; 



+ o(N°) for i > 2 



(54) 



where V2i are the same constant appearing in the case 
of the hard- wall trap [H], i.e. 1/4 = -0.0092552, v 6 = 
0.00404469, etc... 

The above large- N predictions are fully supported by 
the numerical data at fixed N with increasing N, as 
shown in Fig. |9l where we also show data averaged over 
two subsequent particle numbers to suppress the odd- 
even oscillations. In the case of the particle variance, the 
amplitude of the odd-even oscillations appear to decrease 
as 0(1/N), while the average between the data for subse- 
quent particle numbers approaches the predicted asymp- 
totic behavior much more rapidly. In the case of the 
quartic cumulant, the oscillations get suppressed more 
slowly, but their odd-even average approaches the pre- 
dicted value quite rapidly. For the largest available values 
of N the difference of these averages from the asymptotic 
predicted behaviors is O(10~ 5 ). 

We now consider the interval S — [—x, x] around the 
center of the trap. In Fig. [TU] we show the particle vari- 
ance for values of N up to N = 180. They show a be- 
havior analogous to that of the entanglement entropies, 
see Fig. [5j and are consistent with 

Vs (x) « \ [In N + 1 + 7E + In 4 + f v (0] (55) 

The analysis of the data with increasing N is consistent 
with the relation f v (Q = f^ a) (C), where /| q) (C) 

are the 

corresponding scaling functions of the entangelement en- 
tropies, cf. Eq. (|4"gj) . Therefore, fv(C) ls wen approxi- 
mated by the same function f a (0i °f- Eq. (|49)). Again, 
this behavior resembles that of the homogeneous system 
within a hard-wall trap of size L = 21 = 2, which is [l9| 



1 



Vs(x) = —z ln7V + lnsin(7ra;) + 
+l+ lE + hi2 + 0{N~ 1 ) 



(56) 



In Fig. [TT] we show the third and quartic cumulants of 
the interval S — [— x, x]. They are characterized by os- 
cillations which increase when £ — > y/2, but remain ap- 
parently limited with increasing N. 



IV. HIGHER-DIMENSIONAL SYSTEMS 

In this section we consider Fermi gases confined by two- 
and three-dimensional traps. We again study the large- 
TV behavior of the particle correlators, cumulants of the 
particle-number distribution and entanglement entropies 
of extended spatial regions. 
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FIG. 10: (Color online) The particle variance of intervals 
S = [— x,x] for some values of N, versus £ = x/N 1 / 2 . We 
plot tv 2 V — (In N + 1 + 7_b + In 4) The line shows the function 
ED. 
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FIG. 11: (Color online) The third (top) and quartic (bottom) 
cumulants of the interval S = [—a;, x\. 



The vN and Renyi entanglement entropies of ex- 
tended spatial subsystems in the ground state of ho- 
mogenous Fermi gases of d dimension grow asymptoti- 
cally as N( d ~ 1 ')/ d In N, with a prefactor that is analyti- 
cally computed using the Widom conjecture [3!| , for both 
periodic and open boundary conditions. The logarith- 
mic correction to the power-law behavior is related to 
the area- law violation in lattice free fermions [40l - l47j . i.e. 
for a large subsystem A of linear size I in an infinite d- 
dimensional lattice the entanglement entropies scale like 
S {a) {A) - l d ~ x hit In this section we study the effects 
of a space-dependent confining potential in 2D and 3D 
Fermi systems, investigating again the relations between 
particle fluctuations and entanglement entropies. 
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FIG. 12: (Color online) The particle density and the density 
correlator for 2D systems in a harmonic trap: N~ 1 ^ 2 p(r) vs 
r/N 1/4 (top) and N^Gn&x) vs N 1/4 r (bottom) where r = 
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A. Particle density and its correlator 

Using the results of Sec. [Ill we can easily obtain results 
for the particle density and the density correlator in the 
presence of trap. Some data for 2D and 3D systems in a 
harmonic trap are shown in Figs. 1X21 and 1131 They show 
the scaling behavior 

p(x) « N^Rpir/N 1 ^), (57) 
G„(0,x)«7Vi? G (0,r7V 1 /( 2d )), (58) 

where r = |x|. Note that, even in dimensions higher than 
one, the large- TV space rescalings of the particle density 
and its connected correlation are different. 



Half-space entanglement entropies and particle 
fluctuations 



In homogeneous systems with periodic and open 
boundary conditions, the half-space entaglement en- 



tropies of a square L 2 system with open boundary con- 
ditions behave as 11311 



1 



larrVa 



(59) 



The asymptotic large- N behavior of the half-space parti- 
cle cumulants and Renyi entanglement entropies can be 
also computed analytically in the presence of an external 
harmonic potential. For this purpose, we exploit the fact 
that the corresponding overlap matrix ([22]) is a block di- 
agonal matrix. Indeed, relabeling the indeces n, m of the 
N x N overlap matrix as m, rid, using Eq. ([5]), we can 
write the half-space overlap matrix as 



fi ni , . . . ,nrf ;mi ,. . . ,m 



" />oo 

d =TT<Wm, / dz 4> ril (z)4> mi (z) (60) 

i-o JO 



where <j> n are the ID eigenfunctions ([7]), the indeces 
ni,...,n<j correspond to the lowest N states according 
to the Eqs. © and ©. 
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FIG. 14: (Color online) In the top figure we show the half- 
space vN and a = 2 Renyi entanglement entropies of 2D 
systems with harmonic trap. The dashed lines show the 
predicted asymptotic behaviors (|68[) . In the bottom figure 
we plot subtracted entanglement entropies to further check 
the large- N convergence to Eq. (|68[) . which is clearly demon- 
strated by the data. 



Let us first consider a 2D system. We construct the 
ground state of a Fermi gases by filling all states with 



rii + ri2 < n e , rii, n e = 1, 2, 3... 



(61) 



The number N of particles is a function of n e , which 
asymptotically reads N — n 2 /2. Since the overlap matrix 
(f6Tj) is block diagonal, for any integer k we have 



TrA 2D [N(n e )} k 



TrA 1D K-m) fe , 



ni — 1 
o 

Ai D (Af)„ m = / dzcp n (z)<f) m (z), 



(62) 
(63) 



where Ayo(M) is the half-space M x M overlap matrix 
of the ID system. This also implies analogous exact re- 
lations for all observables which can be constructed by 
traces of powers of the overlap matrix or from its eigen- 



values, such as the particle cumulants and the entangle- 
ment entropies, cf. Eqs. flM]) and If2"8"]). Thus, 



S&W^K)] = E 41dK - m), (64) 



ni — 1 



nil 



(65) 



In order to derive their large-iV asymptotic behaviors, 
we replace the sums by integrals and use the relation 
N = 77,2/2, i.e. 



(a) 

HS,2D 



(AO 



/2N 



dnS^(VW-n), (66) 



(m) 



HS,2D 



(N) 



I IN 



(m) 



HS.1D 



(V2N-n). 



(67) 



Then we use the asymptotic formulas for the ID quanti- 
ties, cf. Eqs. (|3T ]) , (|g2" ]l . (|52"J). ([53"f. obtaining 

5 hs, 2 d(A0 ~ aAr1 / 2 [In TV + a + o(N )] , (68) 



1 



a = 



12^2 



a = 2y a - 2 + 7 In 2. 



The approximations used to derive this asymptotic be- 
havior from Eq. (|64l) should not affect the leading 
0{N 1 / 2 InN) and next-to-leading 0{N 1 / 2 ) term, so that 
the constants a and ao should be considered as exact. 
This is confirmed by the analysis of the large- N behavior 
of numerical data at fixed N. In Fig. [14] we compare these 
asymptotic expansions with the data up to N s=s 5000 for 
the vN and a — 2 Renyi entropy, which clearly support 
them. 

For the particle cumulants we obtain 

Vhs,2D (AO = vN 1 ' 2 [In N + v + o{N )] , (69) 
' v = 2j E + 7 In 2, 



2 3 / 2 tt 2 ' 



and 



^HS « V^VrnN 1 / 2 , 777 > 2, 



(70) 



where i/ m are the constants of the leading large- N behav- 
ior in one dimension, cf. Eq. (f54|) . 

The above calculations can be straightforwardly ex- 
tended to higher dimensions. In three dimensions we 
obtain 



S 



(«) 

HS,3D 



(AO 



vi 



iV 2/3 lniV. 



(71) 



The above method can be also used to express the en- 
tanglement entropies and particle cumulants of the sub- 
systems B and S, defined at the end of Sec. Ill B 21 in 
terms of sum of ID contributions for the interval [—00, x] 
and [—x, x] respectively. In particular, in the case of 2D 
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FIG. 15: (Color online) The space dependence of the a = 
2 Renyi (top) and vN (bottom) entanglement entropies of 
2D systems trapped by a harmonic potential. We plot 
N- 1/2 S ( ^(x) vs x/N 1/4 for the harmonic trap. 



0.O 



r < 

J^-0.5 



3D harmonic trap 



v N=56 

° N=165 

* N=455 

x N=1330 

< N=2024 

o N=2600 



0.5 



1.0 



/<5 



FIG. 16: (Color online) The space dependence of vN entan- 
glement entropies of 3D systems trapped by a harmonic po- 
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stripes S contained within two parallel lines at distance 
x from the center, the overlap matrix is 



a ni,...,7i < i;mi, 



i=2 



dz (f> ni (z)(f> mi (z), 



(72) 



which leads to equations analogous to Eqs. (|62H65I) . 
Then, using the continuum approximation and the 
asymptotic large-* behaviors of the ID entanglement 
enetropies of the interval [— x, x], cf. Eq. (|3Hj). we arrive 
at the asymptotic behavior 

s s a L( x ) ~ V2C a N^ 2 \lxiN + h sM (x/Ny 4 )\.(73) 

The coefficient of the leading logarithmic term is just 
twice that of the half-space entanglement entropy ([79]), 
because the boundary of the stripe is double. Analo- 
gously, in three dimensions, we obtain 



7 (°0 

'S,3D 



(,:) 



~7f 



_/V 2 / 3 



lnN + h sia) (x/N^ 6 ) .(74) 



Note that the large- N scaling of the space variables de- 
pends on the spatial dimension d. For a generic d, it 
depends on the scaling variable x/N 1 ^ 2 ^ . 

The large- N scaling of the space dependence of the en- 
tanglement entropies from the size of the extended space 
region can be also checked from the difference Sa(x) = 
Sb (x) — Srs (we recall that the subsystem B is separated 
from the rest by a hyperplane at a distance x from the 
center of the trap). Fig. 1151 shows numerical results for 
the a — 1 vN and the a = 2 Renyi entropyes, which sup- 
port the large-* scaling S { * } (x) = N 1 ' 2 ^ ] (x/* 1 / 4 ), 
i.e. the same large- TV scaling of the space dependence as 
in Eq. (1731 . Fig. [TE1 shows the vN Sa(x) for 3D systems, 
which appear to scale as Sa{x) = * 2 / 3 f&ix/N 1 ^). 

Analogous results are obtained for the particle vari- 
ance. In particular, we find that the ratios of the coef- 
ficients of the leading terms in the asymptotic behaviors 
of the entanglement entropies and the particle variance 
satisfy the universal relation 



Si a) /VA* (1 + a 6 ~ )7r +0(1/ In*) 



(75) 



for any subsystem A considered and in any dimension. 



V. CONCLUSIONS 

We investigate the quantum correlations arising in the 
ground state of free fermion gases trapped by an external 
space-dependent harmonic potential, V cx x 2 jl 2 where I 
is the trap size, in one, two and three dimensions. We 
consider systems of * particles, and focus on the large- 
* scaling behaviors of the quantum correlations, as in- 
ferred by the expectation values of product of local opera- 
tors and bipartite entanglement entropies which quantify 
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the nontrivial entanglement connections between differ- 
ent parts of extended quantum systems. In particular, we 
study the relations between the entanglement entropies 
and the cumulants of the particle distribution within the 
same extended subsystem, which can be obtained by in- 
tegration of the particle-density correlations. 

Our results for the large- iV behaviors of the particle 
density p{x), the two-point particle correlation C{x,y) 
and the connected density-density correlation G n (x,y), 
can be summarized by the following scaling equations: 



p(r) « N e r d R P (N^/ d r/0, r=|x|, (76) 



and 



C(x 1)X2 ) a N e C d Rc{N e ' d *i/ti, JV e/d x 2 /0, (77) 
G„(x!,x 2 ) « N 2e C 2d RG(N e / d xi/Z,N 6 / d x 2 /Z), (78) 

for X! 7^ x 2 , where d is the spatial dimension of the sys- 
tem, £ = l e is the (oscillator) length scale induced by the 
trap, and 9 = 1/2 is the trap exponent for the harmonic 
potential. The above large- N behaviors are expected to 
also hold for higher power laws of the external poten- 
tial, i.e. V(r) oc (r/l) p , by replacing the corresponding 
value of the trap exponent, i.e. 9 = p/(p+2). In the limit 
p — > oo, corresponding to hard- wall trap, the scaling laws 
of homogeneous systems are recovered by setting 9 = 1. 

We compute and analyze the asymptotic large- N be- 
haviors of the particle cumulants and entanglement en- 
tropies of extended spatial regions. Our main results are: 

(i) The half-space Renyi entanglement entropies be- 
have as 

S$ = [In TV + c Q + o(l)} , (79) 

which includes the vN entanglement entropy when a — ¥ 
1. In ID systems, the constant of the logarithmic term is 
equal to that of the homogeneous system, i.e. c; = 1/6, 
which is related to the central charge c = 1 of the cor- 
responding conformal field theory [24], [35| ■ We also de- 
termine the subleading constant cq, cf. Eqs. (|4HI42p . We 
also obtain the constants c; and cq in higher dimensions, 
cf. Eqs. (f6"5f and (I7T|) : in particular we find q = l/(6\/2) 
and ci = 1/ (6-^/6) for the leading logarithmic term in two 
and three dimensions respectively. 

(ii) We compute the asymptotic large- N behavior of 
the half-space particle cumulants. Only even cumulants 
are nonzero, because half-space odd cumulants vanish by 
symmetry. We obtain 

I^hs = v l N^ d - 1 ^ d [In N + v + o(l)] (80) 
Vg k) ^w 2k N^/ d , fc>2, (81) 



In one dimension, see Eqs. (|52"]l and (j5"4")h the constants 
of the leading terms u; and W2k turn out to be equal to 
those of the homogeneous system with open boundary 
conditions (hard walls), which were already computed in 
Refs. [lj| (see also Refs. [ill IHI )• I n particular, vi = 



l/(27r 2 ) and vo is reported in Eq. (l52j) . The constants 
vi and vq are also evaluated in higher dimensions, cf. 
Eqs. ([69]) and (f70j) . Only the particle variance presents 
the leading logarithmic term, like homogeneous systems. 
We find that, in any dimension and for any subsystem 
A, the ratio of the coefficients of the leading terms in the 
entanglement entropies and particle variance satisfies the 
relation 



Vi 



7T 

y 



(82) 



(iii) We also consider spatial bipartitions with different 
geometries, in particular the entanglement entropy of a 
stripe S around the center of the trap with the boundaries 
at a distance x (in ID S — [— x, x]), defined in Sec. lIIB^I 
and studied its space dependence. We find the general 
behavior 



1 



-2ciN {d - 1)/d 



hxN ■ 



(83) 



where cj is the same constant appearing in Eq. (|79p . The 
coefficient of the leading logarithmic term is just twice 
that of the half-space entanglement entropy (f79| , in any 
dimension, essentially because the boundary of the stripe 
is double. A detailed analysis of the ID case is reported 
in Sec. IIIIB 21 The particle variance shows an analogous 
behavior, i.e. 

V w 2viN (d - 1)/d \kxN+ +fv{N (e - 1)/d x/0] , (84) 



where Vi is the same constant appearing in Eq. (|80j) . Note 
that the large- N scaling of the space dependence of the 
entanglement entropies and particle variance is analogous 
to that of the particle density, while it differs from that of 
the particle correlation C(x, y) and the connected density 
correlation G n (x,y), cf. Eqs. ([77)1 and (175]) . 

The above results (i), (ii) and (iii) are consistent with 
the known asymptotic behaviors of homogeneous Fermi 
gases with open boundary conditions [l3|, [l9|, [26| , obtain- 
able by setting = p/(p + 2) 1. 

The large-iV asymptotic behaviors are rapidly ap- 
proached with increasing the number of particles. For 
example, in one dimension the behavior of O(10 2 ), of 
even less, particles is already well approximated by the 
asymptotic behaviors. 

Finally, a few comments are in order concerning the 
relations between particle cumulants and entanglement 
entropies of an extended subsystems A. For noninteract- 
ing fermions, one can write down a formal expansion of 
the entanglement entropies of bipartitions in terms of the 
even cumulants 

EM3. 

such as 



— vf ] + 

945 " 



j(2) 



7i 



-Va 



192 



' A 



23040 A 



(85) 
(86) 
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In homogeneous noninteracting fermion gases with N 
particles in a finite volume of any dimension d, the 
above expansions gets effectively truncated in the large- 
N limit because the high cumulants with m > 2 
are all suppressed relatively to the particle variance Va- 
The leading i\K d_1 )/ d hx./V asymptotic behavior of 



in Eqs. (|85l) and (|86|) arises from Va only, because the 
leading order of each cumulant with k > 2 vanishes 
for any subsystem A (including disjoint ones) in any di- 
mension. This implies the general asymptotic relation 



S 



(Q) (1 + a- 1 ) 7 r i 



(i 



(87) 



Our results for Fermi gases trapped by a harmonic poten- 
tial show an analogous scenario: the asymptotic relation 
(|87f holds as well, in any dimension, significantly extend- 
ing its validity. 

We mention that the close relation between entangle- 
ment entropy and variance in noninteracting Fermi gas is 
also found in off-e quili brium phenomena after local quan- 
tum quenches (l5l Il8l |49| , and in some dynamics regime 
of the off-equilibrium expansion of Fermi gases from a 
trap [50j . 

The situation is more involved for interacting systems. 
In systems with localized interactions arising from im- 
purities [5l|, all the cumulants y( 2fc ) contribute to the 
asymptotic large- N behavior of the entanglement en- 
tropies in the expansion 



(a) 



Eoo 
fc=l 



al- 



though the expansion turns out to be rapidly converg- 
ing [1{|. The conservation of a global charge, and in 
particular the particle number, is crucial for the connec- 
tions between bipartite entanglement and particle fluctu- 
ations. For interacting systems not conserving the par- 
ticle number, the entanglement should be related to the 
more fundamental energy transport. 
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Appendix A: Asymptotic behavior of the ID 
half-space entanglement entropies 

We consider a ID lattice model of spinless fermions in 
the presence of an external space-dependent potential 



where Cj is a spinless fermion operator, and 



(Al) 



h l3 = J (Sij - - 7^J+i) + V(xi)5ij. (A2) 



We consider a power-law spatial dependence for the trap- 
ping potential, 



V(r) = -v p r p , 
P 



(A3) 



where r is the distance from the center of the trap, v is 
a positive constant and p an even integer number. The 
trap size I is defined as I = J 1 / p /w. In the following we 
set J = 1 . 

In one dimension this lattice free-fermion model can 
be exactly mapped into the hard-core (HC) limit of the 
Bose-Hubbard (BH) model, see, e.g., Ref. [521 ] . 



J 



(A4) 



(ij) 



where (ij) is the set of nearest-neighbor sites, hi are 
bosonic operators, rii = b\bi is the particle density op- 
erator, and TV = Q^. n,*) is the particle number. The HC 
limit U — > oo of the BH model implies that the particle 
number rij per site is restricted to the values n$ = 0,1. 
As a consequence of their exact mapping, the HC BH 
and lattice free spinless fermions share the same particle 
density and density-density correlation, particle distri- 
bution of extended space regions, and also entanglement 
entropies of connected spatial bipartitions (l^ |. 

The large-Z limit keeping N fixed differs from that per- 
formed at fixed chemical potential /i, i.e., considering the 
BH Hamiltonian 



Hp — "Hbh 



(m-i)E ; 



(A5) 



Indeed, the large trap-size limit, keeping /i fixed, implies 
an increase of the particle number so that 



N/l d = fa) 



(A6) 



asymptotically, where d is the spatial dimension and p(/i) 
is a finite function of /i. This thermodynamic limit is 
usually considered when quantum transitions are studied 
in confined particle systems. In the absence of the trap, 
the ID HC BH model has three phases: the empty state 
for > 1 with (rii) = 0, which may be seen as a particular 
n = Mott phase, a gapless superfluid phase for \fi\ < 1, 



and a n = 1 Mott phase for fi < — 1. See, e.g., Ref. [52[. 

We consider chains with even sites L and open bound- 
ary conditions, and a trap of size I centered between the 
middle sites of the chain. We divide the chain in two 
connected parts of length I a and L — I a, and consider 
their Renyi entropies 



S (a \lA;L) = S^(L-l A ;L) 



1 



1 - a 



lnTrp A (A7) 



where pa is the reduced density matrix of one of the two 
subsystems. Let us consider the half-space entanglement 



(a) 
HS 



5 (Q) (L/2,L) 



(A8) 



14 



Its large-L behavior can be written as 

\nL + e a + 0{L- 1/a ) 
l + a- 1 



c(°0 _ r< 
°hs — °<> 



12 



2 Q = In y/l - fi 2 + ln(4/7r) + y a , 



Va 



50 dt r 

T 
l 

sinht 



1 - cr 



asinhi/a; sinht 



(A9) 

(A10) 
(All) 
I x 

(A12) 



where y a is given in Eq. (|A12I) . This equation has been 
obtained using the results of Refs. [22|, |35|, l53l. |54|. 
In a system confined by a trap of size I, we have [33| 



?(<*) = 
>hs — 



C a 



Um L ^ 00 S^(L/2;L) = (A13) 



\nCe+e a + 0(C 1/a ) 



where 



£ e = a e (n)l 



(A14) 



is the entanglement length, which also enters the asymp- 
totic formula of the energy difference of the two lowest 
states 



A 



t{4>), i(0) = 1/2 -1^-1/21, (A15) 



where <j) is phase-like variable, < <f> < 1, which 
parametrizes the modulations of the amplitude due to 
the periodic asymptotic occurrence of level crossings in 
the large-? limit |36j ] . 

We now consider a HC BH system of N particles in a 
trap (centered in the middle of the chain), and we want 
to study the bipartite entanglement entropy, and in par- 
ticular the half-lattice entanglement entropy, in the limit 
of large trap size (after taking the infinite chain limit 
L — > oo ) as a function of N. More precisely, we consider 

S^(l,N)=Um L ^ 00 S^(L/2,L;l,N), (A16) 



'HS 



(A17) 



which is a finite function of N. Since N is kept fixed while 
I — > oo, this corresponds to the dilute regime [i < 1. Let 
us define S = 1 — fi, thus S — > + corresponds to the 
dilution limit. 

We want to derive asymptotic large-TV behavior of 
Shs(N). For this purpose, we need: (i) the dependence 
of a e on [i when /i — > 1~; (ii) the relation between the 
number of particles and the chemical potential. 

The dependence of a e on fi can be inferred from the 
behavior of the gap for 5 — > (see Sec. IV B of Ref. 3(f]), 
by matching Eq. (| A15|) with its asymptotic behavior for 



(A18) 



where 



c p (5) 



2V7Fr(3/4) 

r(i/4) 

7T 

V2 



for p = 2 
for p = 4 
for p — > oo 



(A19) 



In ID particle systems, the thermodynamic limit at 
fixed fx corresponds to N, I oo keeping the ratio AT/Z 
fixed. Indeed, we have 



N 



(A20) 



The function p(/i) can be computed in the HC limit. The 
particle density in the large-/ limit turns out to approach 
its local density approximation (LDA), with corrections 
that are suppressed by powers of the trap size and present 
a nontrivial TSS behaviour. Within the LDA, the parti- 
cle density at the spatial coordinate x equals the particle 
density of the homogeneous system at the effective chem- 
ical potential 



1 fx \p 
H cS (x) = fi+- [j) 



(A21) 



The LDA of the particle density reads (n x )ida = 
p\A<x{x/l), where 

!0 for n e s(x) > 1, 

(1/n) arccos^ cff (x) for -1 < n e s(x) < 1, 
1 for fi cS (x) < -1. 

(A22) 

Asymptotically, the total particle number is obtained by 
integrating the LDA of the particle density pida, obtain- 
ing 



Pida(y) dy. 



In the low-density regime, S = 1 — jj, — > 0, 
p(/x) = r p ^+ 2 )/^[l + 0(5)], 



with 



1 for p = 2 

EWf) f 0r D = 4 

3V^r(3/4) 101 P * 

for p — > oo 



(A23) 



(A24) 



(A25) 



Using these results we obtain the entanglement entropy 
(|A13I) in terms of N: 



S$(N)=C a 



In N + In 



(A26) 



The above formula can be compared with large-L and 
then large-/ extrapolations of numerical results by exact 
diagonalization of chains of size L and traps of size I 
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centered between the middle sites. In particular, the TV- 
dependent limit (|A26[) is approached with 0{l~ 26 ) cor- 
rections, where 9 — p/(p + 2) is the trap exponent, con- 
forming gen eral theoretical arguments on the large- 1 cor- 
rections [32i | . 

The large trap-size limit keeping the particle number N 



fixed corresponds to the dilute limit of continuum models 
in the presence of the trapping potential. As argued in 
Ref. (321, it also represents the asymptotic large- N be- 
havior for finite on-site couplings U > within the BH 
model, or finite-strength models of bosonic gases with 
short-range interactions. 
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